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Abstract 

We develop a recipe which allows one to recursively and uniquely decouple the 
large scale from the small scale in mass matrices of the seesaw type, up to any 
order in the inverse of the large scale. Our method allows one to calculate the mass 
matrix of the light neutrinos with arbitrary precision. The same method can be 
applied in the case of quark mass matrices in an extension of the Standard Model 
with vector-like quarks which have mass terms at a scale much higher than the 
electroweak scale. 



1 Introduction 

Neutrino physics might prove to be the first window to physics beyond the Standard 
Model. The reason for this is the evidence for neutrino oscillations found, in particu- 
lar, in experiments measuring the atmospheric neutrino flux Most probably, physics 
beyond the Standard Model is associated with a new mass or energy scale above the 
electroweak scale. This idea has been exploited in the seesaw mechanism to explain the 
smallness of the neutrino masses relative to the charged-lepton masses 0. Inspired by 
the indication for non-zero neutrino masses from neutrino oscillations, model building 
for incorporating neutrino masses and lepton mixing and for explaining features of the 
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neutrino mass spectrum and of the mixing matrix is in full swing. In this context, the 
seesaw mechanism plays an important role 0]. 

Usually, the lowest order in the inverse heavy seesaw scale is sufficient for the purpose 
of describing the light-neutrino mass matrix. However, one can envisage situations which 
warrant going beyond the lowest order, like 

• the need to know the neutrino masses to a better precision; 

• a degeneracy of the mass spectrum at lowest order, which is lifted at higher orders; 

• the presence of massless neutrinos at lowest order, which become massive at higher 
orders; 

• a lepton mixing matrix with some features at lowest order, which disappear at higher 
orders 0; 

• light mass terms which, though much smaller than the heavy ones, are not all of 
the same order of magnitude; or, inversely, heavy masses with different orders of 
magnitude ||^; or the same situation for both types of mass terms. 

In this letter we propose a simple, minimal procedure which allows one to decouple the 
small seesaw scale from the large one order by order in the inverse large scale. The 
procedure is defined by a particular Ansatz for the unitary matrix W which performs 
this decoupling, together with a series expansion for W in the inverse large scale. This 
procedure enables, in principle, the calculation of W, of the mass matrix of the light 
neutrinos, and of the mass matrix of the heavy neutrinos, up to any arbitrarily large 
order in the inverse heavy mass scale. 

Seesaw-type mass matrices, with two different mass scales, are not confined to neutrino 
physics. In this letter we also show how one may recursively decouple the light from the 
heavy quarks when one adds vector-like quarks to the Standard Model and one assumes 
the gauge-invariant mass terms to be very large. 

In Section |^ we develop the decoupling procedure in the case of the usual seesaw 
mechanism for neutrinos, whereas in Section |^ we consider the extension of the Standard 
Model with vector-like quarks. In Section ^ we present our conclusions. 

2 Neutrinos 

Let us consider a general model in which, to the generations of leptons of the Standard 
Model, one adds right-handed (singlet) neutrinos z//j One may take the charge con- 
jugates of the latter fields, z/^ = CVr^ , which are left-handed, and put them together with 
the doublet left-handed neutrinos. (The matrix C is the Dirac-Pauli charge-conjugation 
matrix.) The most general mass terms for all the neutrinos fields are then 




The ul X ul matrix and the nn x nji matrix Mr are symmetric; the matrix is 
X riL and arbitrary. 
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The mass terms in M^) are | AJ| = 1/2, where / denotes the weak isospin, and originate, 
usually, in the Yukawa couplings of the leptons to scalar doublets which acquire a vacuum 
expectation value (VEV). The mass terms in are |A/| = 1 and are present if there 
exists in the theory either a fundamental or an effective Higgs triplet with a VEV. The 
mass terms in Mr are gauge-invariant and we shall assume them to be much larger than 
the ones in either Mo or Ml. We shall moreover assume Mr to be non-singular. The 



order of magnitude of the eigenvalues of y Mj^Mji will be denoted by rriR. We shall 
assume that all matrix elements of Ml and of Md are much smaller than rriR. Under 
these assumptions, there are in the theory ur heavy neutrinos with mass of order tur, 
and hl light neutrinos with masses which are, in the case Ml = 0, suppressed by one 
or more inverse powers of rriR |]^.[| Note, however, that the notion "order of magnitude 
tur" should not be taken too literally. To allow for a more general situation it should 



be understood in the sense that the eigenvalues of y Mj^Mr, though not necessarily all 
of the same order of magnitude, are all much larger than the matrix elements of Ml and 
Ml)] the latter, too, may have orders of magnitude different among themselves. 

It is our aim to decouple the heavy from the light neutrino fields, and to derive 
the effective mass matrices for each of them. This is done by performing a unitary 
transformation of the neutrino fields by means of an {ul + ur) x [ul + ur) unitary matrix 



W\ ^"^^^ ) , (2) 



\ ) \ ^heavy 

and by requiring from this unitary transformation that 

where Mught and Mheavy are, respectively, ul x ul and ur x ur symmetric matrices. 

In Eq. (|^) we require that the transformation performed by the matrix W has the 
property of leading to a zero ul x nR submatrix of the mass matrix. On these gen- 
eral grounds, W must have nLUR degrees of freedom, and any degrees of freedom in W 
beyond these ones are not needed and might even render its computation unnecessarily 
complicated. We shall therefore make the following y4nsat20 for W: 

/ Vl - BB^ B \ 

where B is an ul x tir matrix which must be fixed as a function of Ml, Md, and Mr. 
The square roots in Eq. (H) should be understood as power series, for instance. 



1 _ BB^ = 1 - iBB^ - \BB^BB^ - ^BB^BB^BB^ (5) 



"singular seesaw mechanism" [Q, where the matrix Mr is singular, is also possible. If \J m\^Mr 
has rto eigenvalues 0, then we can define n'j^ ^ ul + ng and rt^ ^ nR — hq, and this singular seesaw case 
in contained mathematically in the usual seesaw case with n'j^ left-handed doublets and n'^ right-handed 
singlets. This is most easily seen by going into a basis where Mr is diagonal, with the uq diagonal zero 
entries coming first. 

^In Ref. jg] an equivalent Ansatz has been proposed; Ref. ^ uses the Ansatz of Eq. in a more 
restrictive context. 



3 



The matrix W in Eq. (^) is unitary by construction. Indeed, Eq. (Q) might be seen as a 
generahzation, for matrices, of the usual form of a 2 x 2 orthogonal matrix 

/ Vl — sin^ 6 sin 6 \ 



sm U V i — sm 



With the Ansatz for W in Eq. (H), the condition of the vanishing of the off-diagonal 
submatrices in Eq. (^) reads 



B'^Ml\Ji - BB^ + VI - BTB*Md\/1 - BB^ - B^MJ^B^ - sjl - BTB*MrB^ = 0. (6) 
This equation may be solved by assuming that 5 is a power series in l/mn, namely, 

B = B1 + B2 + B3 + B4 + ..., (7) 



1-BB^ = 1 - ^BiBl - ^ [BiBl + B2BI 

-i {b,bI + B2BI + B,Bl + \B^BlB^Bl) (8) 

where Bj is by definition proportional to {mn)'-' . The recursive solubility of Eq. can 
be seen in the following way. At order {niRf Eq. simply reads Md — MrBI = 0, and 
this fixes Bi . We then note that the i-th. order of Vl — BB^ is a function of the Bj with 
j < i, cf. Eq. (^. Following this observation, an inspection of the k-th order of Eq. (^) 
shows that, except for a term —MfiBl_^^, all terms only involve the Bj with j < k. In 
other words, the k-th order terms of Eq. lead to an expression for -B^+i in terms of 
the Bj with j < k. This proves the recursive calculability of B. 
From Eq. (|^) one thus obtains 

MrBI = Md, (9) 
MnBl = M],'*M*oMl, (10) 
MrBI = M]^^*M^^MdMIMl- M^^*M*jjMIM^^Md-\MdM\jM^^*M^^Md. (11) 

The expressions for the B^ with k > A are rather long and complicated, but no funda- 
mental problem arises in their recursive computation. The effective mass matrices for the 
light and heavy neutrinos are then 



-MIM^^M^^*M^^MdMIMl - MlMIM]^M^^M^^*M^^Md] + . . . , (12) 
+ i (MdMIMIM^^Mr^* + M^^*M^^MdMIMI) + . . . , (13) 



respectively. Equations ([I^ ) and ( pTSf ) are correct up to order (m/j) ^ and (m^j) ^, respec- 
tively. 
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In the most common case of the seesaw mechanism, which is for instance obtained by 
extending the Standard Model with right-handed neutrino singlets but avoiding the pres- 
ence of any Higgs triplet, one has a vanishing matrix M^. Then the recursive decoupling 
mechanism becomes much simpler because 

Mi = ^ B2 = Bi = Be = ... = 0. (14) 



Let us sketch a proof of Eq. (|14[) by induction. By direct calculation we find that B2 = 
when Ml = 0, see Eq. (|10D. Now we assume that B2 = = . . . = = {i > 1). In 
order to calculate B21+2 from Eq. (^, we need to know all the Bj with j < 2£ + 1. Then, 
according to the induction assumption, only the Bj with odd index j contribute to -82^+2- 
We now notice that Vl — BB^ only has even powers of l/mn up to order 2i + 1. Then, 
a term-by-term inspection of Eq. (P) with Ml = reveals that all terms are zero at order 
2i+l except for the term — M^i^l^^g- Therefore, B21+2 = 0, and this proves the theorem 
in Eq. As a consequence, we have 

B = B^ + B, + ..., (15) 



l-BB^ = l-lBiBl-UBiBl + BsBl + lBiBlBiBl)-..., (16) 



with 

MrBI = Md, (17) 
MrBI = -M^^*M*j^MIM^^Md-\MdM\jM^^*M^^Md, (18) 

and, using similar considerations as in the above proof by induction, we also find that 
only odd powers of 1/m/j are present in M^ght and Mheavy, in particular, 

Miight = -MIM],^Md 

+\MlM],^ (MdMIM],^* + M^'*M*dMI) M^'Md + ..., (19) 

Mheavy = Mr + ^ {MdMIM^^* + M]^'* M*r,Mj,) + . . . (20) 

It is possible to proceed and derive, recursively, expressions for B^, B^, and so on, and 
also better approximations for Mught and for Mheavy However, that task is tedious and 
complicated, since the number of terms in the polynomials increases very rapidly. 

The subcase of the case Ml = where the leading term in Mught (|19|) is zero has been 
studied in Ref. 0. In this subcase one does not need recursive expansions, as one is able 
to give closed forms for all the relevant matrices. By using Eq. (^ and the definition 
Z = M]^^Ml,, one readily verifies that P 

Mi = , MIM^^Md = ^ B = zUl + ZZ^) '^^^ = (l + Z^Z) '^^^ Z^ . (21) 



Then, using Eq. (|^), one finds after some algebra that 

Mught = and Mheavy = (l + Z* Z'^) Mr (l + ZZ^) . (22) 

It is remarkable that the vanishing of the lowest order contribution to Mnght in Eq. (|T9|) 
entails the exact vanishing of the whole matrix.| 

■^We thank A. Pilaftsis and A. Joshipura for calling our attention to this interesting subcase. 
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3 Vector-like quarks 



Let us consider that to the Ug generations of u and d quarks of the Standard Model one 
adds n' vector-hke singlet quarks of, say, charge —1/3. One then has left-handed quarks 
(11, which are components of doublets of SU(2), together with n' left-handed quarks Dl, 
which are singlets of SU(2); the Ug + n' right-handed quarks dji are all singlets |]10[- The 
mass term for these down-type quarks reads 



4, DL)MdR + H.c, with M={ „ . (23) 




The matrix m is x {jig + n') and gives the |A/| = 1/2 mass terms; the matrix M is 
n' X {jig + n') and contains the SU(2)-invariant mass terms. We shall assume that the 
matrix elements of m are of order A, and that the n' eigenvalues of yMM^ are all non-zero 
and of order A, with A/A = e <C 1. 

The squared- mass matrix for the left-handed quarks d^ and is 

« .< « ( rnm) mM^ \ 

= [ Mm^ MM^ ) ■ (24) 

After diagonalization, there will be "light" quarks with masses of order A, and n' 
"heavy" quarks with masses of order A. One would like to decouple the light quarks from 
the heavy ones and to find separate squared-mass matrices for each category. This we 
shall do by means of a unitary transformation. 



U^MM^U = '^^;f^* ^ ^ , (25) 




where U is {ug + n') x {rig + n') unitary, while -Dught and -Dheavy are, respectively, Ug x Ug 
and n' x n' Hermitian (and positive definite) matrices. 

In Eq. (^) we are, as a matter of fact, transforming an ng x n' submatrix of A^A^^^ 
to zero by means of a unitary transformation. It is therefore convenient to choose for 
the unitary matrix U of that transformation an Ansatz which has exactly UgTi' (complex) 
degrees of freedom. We therefore write 

(26) 

where F is an arbitrary Ug x n' matrix, to be fixed as a function of m and of M. The 
square roots in Eq. ( p6D should be understood as power series, just as in Eq. (|^). It is 
then obvious that the matrix U in Eq. ( PBj ) is indeed unitary; it is the analogue of the 
matrix W in the previous section. 

The condition which fixes F is the demand that there is a zero submatrix in Eq. (pSf) , 
namely. 




1 - FF^mm^F + Jl - FF^mM^Jl - F^F-FMm^F-FMM^Jl - F^F = 0. (27) 
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We solve this equation by writing F as a power series in e, 



F = F^ + Fs + F5 + ..., (28) 

where Fj has matrix elements of order . Notice that only terms with odd j are necessary. 
Defining 

S=[MM^y^ = Sl (29) 

it is easy to find 

Fi = mM^S, (30) 
F3 = m (m^mM^S - M^SMm)mM^S - S'^Mm)mM^) S. (31) 

Unfortunately, the expressions for the Fj become quite complicated when j > 5, although 
they may in principle be derived recursively without any problem. 

One is thus able to decouple the light from the heavy quarks. One finds their respective 
squared-mass matrices to be given by 

Aight = m(l- M^SM^m^ 

-\m [(1 - M^SM) w)mM^S'^M + M^S'^Mm^m (l - M^SM)] m) 
+ ... (32) 
Dheavy = M + \ [SMw)mM^ + Mm^mM^S) + ... (33) 

The expressions of -Dught and -Dheavy have been given to sub leading order in e^. Note that 
only even powers of e appear in the expansions of these squared-mass matrices. 

It is interesting to write down the equations above in a specific simple weak basis. 
By performing a unitary mixing of the right-handed singlet quarks (Ir and also of the 
left-handed singlet quarks -Dl, one finds a weak basis in which 



G, J ) and M = ( 0, M ) , (34) 



where M is a diagonal and positive n' x n' matrix, while G and J are general Ug x ng and 
Hg X n' matrices, respectively. In this weak basis one finds 

Aight = GG^ -\[GG\JM-^J^] + (35) 

Dheavy = + \ (M'^JUM + M jU M~^) + . . . , (36) 

where {X, Y} = XY + YX is the anticommutator. In this basis it is also obvious that 
the {ug + n') X (n^ -|- n') matrix M^S'M, which appears in Eq. (^), is a projector, which 
projects into the space characterized by the last n' indices. Therefore, at lowest order, 
in the matrix -Dught of Eq. (^) only the left Hg x Ug submatrix of m contributes. The 
degrees of freedom beyond the Standard Model, which are characterized in the weak basis 
of Eq. (Rl) by the last n' indices, do not contribute to the light-quark masses, as expected. 
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4 Conclusions 



In this note we have discussed a method which allows to separate, order by order in 
the inverse seesaw scale, the small and large scales in seesaw-type mass matrices, and 
to obtain a recursive expansion for the light-neutrino and for the heavy-neutrino mass 
matrices. This has been achieved by the Ansatz of Eq. (^) for the matrix W, which 
approaches the unit matrix when the large scale tends to infinity. A similar method, with 
the matrix U of Eq. (p6D, can be used to decouple the light quarks from additional heavy 
quarks in extensions of the Standard Model with vector-like quarks. The lepton mixing 
matrix and the quark mixing matrix can then be calculated from the matrices Mught of 
Eq. (|T^) and Dugtit of Eq. (|3^) , respectively. For this purpose, the lowest order of these 
matrices might be sufficient, yet in some instances it might be necessary to include further 
terms in order to obtain trustworthy representations of the mixing matrices. If one wants 
to go beyond lowest-order mixing, then in the neutrino case there are contributions from 
higher orders in Mught and in W; the same holds in the vector-like-quark case for -Dught 
and U. 
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